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Abstract 

We report on a measurement of the charmless semileptonic B decays, — > 7r~£"'"z/ and B^ 
p''t^u, based on 140 fb~^ data collected with the Belle detector at the KEKB e'^e~ asymmetric 
collider. In this analysis, accompanying B mesons are reconstructed from semileptonic B D^*^iv 
decays, which enables us to detect the signal with high purity and with marginal statistics. We 
found branching fractions of B{B^ Tr-i+u) = (1.76 ± 0.28 ± 0.20 ± 0.03) x lO"'' and 6(5° 
p~l'^u) = (2.54ib0.78±0.85±0.30) x 10~^, where the errors are statistical, experimental systematic, 
systematic due to form-factor uncertainties. We present also the branching fractions in three 
intervals; < 8, 8 < < 16, > 16GeV^/c^, for each decay mode. Based on these results, 
the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element Vub is extracted. All of the 
presented results are preliminary. 

PACS numbers: 12.15.Hh, 12.38.Gc, 13.25.Hw 
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INTRODUCTION 



Exclusive B Xud-v decays proceed via h — > uW~ tree diagram with a spectator quark, 
and can be used to determine \Vub\i one of the smallest and least known elements of the 
Cabibbo-Kobayashi-Maskawa matrix The major theoretical issue is determination of 
the form-factors (FF) involved in the decays. 

In principle lattice QCD(LQCD) provides the most reliable calculation of FFs, but so far 
only quenched calculations have been available. Recently preliminary results of unquenched 
calculations have become reported 0, |3| so that a model-independent determination of 
\Vub\ becomes feasible. Since LQCD results are available only in the high region the 
clean extraction of the B -niv branching fraction in the high (f region (> 16GeV^/c^) is 
essential. 

There have been several measurements in the past by CLEO, BaBar and Belle for the 
B — > vrfz/, p£z/, r]iv and ujiv modes 0, IE B Si- ^ these analyses utilize the method, 
originally developed by CLEO, where the B decays are reconstructed by inferring the un- 
detected neutrino momentum from the missing energy and momentum ("z/-reconstruction 
method") p:]. In the 5-factory era, we may quickly improve the statistical precision by sim- 
ply applying the i/-reconstruction method. However, the systematic uncertainty may soon 
limit the experimental uncertainty because of the poor signal-to-noise ratio. 

In this paper, we present measurements of 5° — > vr'f+z/ and B^ p~£+z/ decays by 
D^*^£u decay tagging, where we reconstruct the whole decay chain from the T(4S'); T(4S') 
BsigBtag, B gig n~{p~)i^u aud Btag ^ D''*^^ i~ u with several sub decay modes. 

The back-to-back correlation of the two B mesons in the T(45') rest frame allows us to 
constrain the kinematics of the double semileptonic decay. Belle has presented a preliminary 
measurement of the inclusive B — > X„£z/ decay using this method [9|. While this method 
gives a signal detection efficiency that is significantly lower than the i/-reconstruction method 
because of the tagging, it allows us to detect the signal with high purity, even with a relaxed 
cut on the lepton momentum down to 0.8 GeV/c. This leads to less experimental systematic 
uncertainty and promises to yield the best overall precision from the large data sample being 
accumulated by 5-factories over the next few years. We present the simultaneous extraction 
of B^ -n-^i^u and B p^i^u and their distribution. Inclusion of the charge conjugate 
decays is implied throughout this paper. 



DATA SET AND EXPERIMENT 

The analysis is based on data recorded with the Belle detector at the asymmetric e^e~ 
collider KEKB operating at the center-of-mass (cm.) energy of the T(4S') resonance jlol |. 
The T(4S') data set used for this study corresponds to an integrated luminosity of 140 fb~^ 
and contains 152 x 10^ BB events. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50- layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals 
(ECL) located inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. 
An iron fiux-return located outside of the coil is instrumented to detect mesons and to 
identify muons (KLM). The detector is described in detail elsewhere jlll |. 

A detailed Monte Carlo (MC) simulation, which fully describes the detector geometry and 
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response and is based on GEANT is applied to estimate the signal detection efficiency 



and to study the background. To examine the FF dependence, MC samples for the 5° 
7r~(p~)i~^i> signal decays are generated with different form-factor models; a quark model 
(ISGW II 111), light cone sum rules (LCSR; Ball'Ol for 7r-£+z/ Q and Bair98 for p-i+u 



15 1) and a quenched lattice QCD (UKQCD |16|). A relativistic quark model (Melikhov 
17[) is also used for p~i~^i>. To model the cross-feed from other B X^d-v decays, MC 



samples are generated with the ISGW II model for the resonant components and the DeFazio- 
Neubert model for non-resonant component {ttIv and plv components are excluded in 
this sample). A MC sample for the D*^l^v decay, to be used for calibration as 

described later, is also prepared. To model the BB and continuum backgrounds, large 
generic BB (based on QQ98 ^) and qq Monte Carlo samples are used. 



EVENT RECONSTRUCTION AND SELECTION 

Charged particle tracks are reconstructed from hits in the SVD and CDC. They are re- 
quired to satisfy track quality cuts based on their impact parameters relative to the measured 
profile of the interaction point (IP profile) of the two beams. Charged kaons are identified 
by combining information on ionization loss (dE/dx) in the CDC, Cherenkov light yields 
in the ACC and time-of-flight measured by the TOF system. For the nominal requirement, 
the kaon identification efficiency is approximately 88% and the rate for misidentification of 
pions as kaons is about 8%. Hadron tracks that are not identified as kaons are treated as 
pions. Tracks satisfying the lepton identification criteria, as described later, are removed 
from consideration. 

Neutral pions are reconstructed using 7 pairs with an invariant mass between 117 and 
150MeV/c^. Each 7 is required to have a minimum energy deposit of E^^ > 50MeV. Kg 
mesons are reconstructed using pairs of charged tracks that have an invariant mass within 
±7.6MeV/c^ of the known Kg mass. 

Electron identification is based on a combination of dE/dx in CDC, the response of ACC, 
shower shape in ECL and the ratio of energy deposit in ECL to the momentum measured by 
the tracking system. Muon identification by KLM is performed by resistive plate counters 
interleaved in the iron yoke. The lepton identification efficiencies are estimated to be about 
90% for both electrons and muons in the momentum region above 1.2 GeV/c. where leptons 
from the prompt B decays dominate. The hadron misidentification rate is measured using 
reconstructed Kg tt+tt^ and found to be less than 0.2% for electrons and 1.5% for muons 
in the same momentum region. 

For the reconstruction of Btag — > D^*^'^i~i', the lepton candidate is required to have 
the right sign charge with respect to the D meson flavor and the laboratory momentum 
greater than l.OGeV/c (p^"^ > 1.0 GeV/c). The candidates are reconstructed by 

using four decay modes of - K^ti^ti^, K^ti^ti^tt^ , KgTT~^ and i^^vr^Tr*^ - and 

six decay modes of - K^tt^, K^iT~^7r^, K^tt+tt+tt^, Kgir'^, Kgir^ir^, Kgir^TT^TC^ . The 
candidates are required to have an invariant mass rriD within ±2a {a is a standard deviation) 
of the nominal D mass, where the mass resolution a is dependent on the decay mode. D*~^ 
mesons are reconstructed by combining the D° candidate and a pion, D*^ D^n^ and 
D*~^ D^Ti^ . They are required to have the mass difference Am = m/)^ — m^, within ±2(j 
of the nominal values. 

For the reconstruction of Bsig — ^ X~d^u, the lepton candidate is required to have the 
right sign charge with respect to the system and p^^^ > 0.8 GeV/c. The X^ system 
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may consist of either a single charged pion or a charged pion and neutral pion candidate, 
and the event may contain no additional charged tracks or 7r° candidates {N^^- = 1 and 
NttO < 1). We also require that the residual energy from neutral clusters is less than 0.3 
GeV {Eneut < 0.3 GeV). The two leptons on the tag and the signal sides are required to 
have opposite charge. Loss of the signal due to — mixing is estimated by the MC 
simulation. 

We then impose a constraint from kinematics of the double semileptonic decay in the 
T(45') rest frame. In the semileptonic decay on each side, -Bi(2) — ^ ^i(2)i^i(2) (^i = D^*'>^i~ 
and Y2 = the angle between the -81(2) meson and the detected Yi(2) system 0B^2) 

is calculated from the relation, pi = {ps — PyY = and the known Pb (the absolute 
momentum of the mother B meson). This means that the -Bi(2) direction is constrained 
on the surface of a cone defined with the angle Ob^^) around the direction of the Yi(2) 
system, as shown graphically in Figure ^ Then the back-to-back relation of the two B 
meson directions implies that the real B direction is on the cross lines of the two cones 
when one of the B system is spatially inverted. Denoting 612 the angle between the D*^i~ 
and the X^i^ systems, the B directional vector ub = {xbjVb, zb) is given by, zb = cos6bi, 
Ub = {cos9b2 — cos6'B2Cos6'i2)/sin6'i2, and 



xb = ±1/1 ;^(cos26'b, + co?P'9b2 — 2cos6b,cos6b2Cos6i2) (1) 

V sm6'i2 

with the coordinate definition in Figure H If the hypothesis of the double semileptonic 
decay is correct and all the decay products are detected except for the two neutrinos, x% 
must range from to 1. Events passing a rather loose cut > —2.0 are used for the signal 
extraction in the later stage of the analysis. 




FIG. 1: Kinematics of the double semileptonic decay. 



In this analysis, Eq. ^has two solutions and the direction of the B meson is not uniquely 
determined. We calculate, therefore, as = (-Efeeam ~ ^Xu)"^ ~ P*xJ'^ using the beam 
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energy {E^^^^), energy {E*x^) and momentum {px^) of the X„ system and neglecting the 
momentum of the B meson in the cm. system. The signal Monte Carlo simulation predicts 
that the average resolution is approximately 0.75GeV^/c^. 

According to Monte Carlo simulations, the major backgrounds come from B^B^ events 
and other B^ — > X^C^v decays, where some particles escape detection. Other backgrounds 
from B^B~ events, B^ X^i'^u decays and qq processes are small. 

With the event selection described above, the signal MC simulation indicates that the 
total detection efficiency (etotai), on the average of the electron and muon channels, is 1.5 x 
10-3 for n~i+u and 8.1 x 10"^ for p'i+u with the LCSR model (Ball'Ol for 7r-£+z/ and Bair98 
for p~£+z/). Here, etotai is defined with respect to the number of produced B'^B^ pairs, where 
one B decays into the signal mode, and includes loss of the signal due to B^ — B^ mixing. 
Because of the relaxed lepton momentum cut (> 0.8GeV/c), the variation of the efficiency 
with different FF models is relatively small; the maximum deviation from LCSR is —1.4% 
with ISGW II for n-i+u and +5.5% with UKQCD for p'i+u. 

The validity of the method to reconstruct the double semileptonic decay is checked by 
using the decay, B^-^ D*~£+z/ followed by D*~ D^tt~,D'^ K~^7r~, with the same 
requirement on the tagging side. Figure |21-a) shows the obtained Mk-ktt distribution and its 
comparison to the MC expectation. With the 140 fb~^ data sample, 147 ± 12 decays are 
clearly identified, while 165 ±9 events are expected based on the product branching fraction 
B{B^ D*-i+u,D*- DH-,DO K+n-) = (1.40 ± 0.07) x lO'^ deduced from 0. 
Their ratio R = 0.89 ± 0.08 is used to correct the total detection efficiency etotai for the 
7r(p)£z/ signals predicted by the signal MC simulations. Figure |21-b) shows the comparison 
of the reconstructed distribution using the same sample to the MC simulation. The 
agreement between the data and MC demonstrates the validity of the present measurement. 




FIG. 2: Reconstructed M{KTnr) distribution(left) and distribution(right) for the D* i^v 
calibration decay. Open circle is data and histogram is signal MC. 
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EXTRACTION OF BRANCHING FRACTIONS 

The — i> Tc~i~^i' and p^i^u signals are extracted by fitting the obtained two- 

dimensional distribution in (x^ , Mx ) , where Mx is the invariant mass of the Xu system. 
Here, the fit components are the two signal modes; B^ 7r^£+z//p~£"^i^, the other B^ — *■ 
X~£~^i' background and the BB background. The PDF (probability distribution function) 
for each component is deduced from the MC simulation. The t^' {p~)l'^v signal events exhibit 
characteristic behavior in both of their x\ and Mx distributions; other B^ — ^ X~^'^v events 
exhibit a weak peaking structure in x\ but a relatively flat distribution in Mx'-, the BB 
background events are distributed uniformly in both variables. The total inclusive branching 
fraction B{B — > is constrained to be 0.25% j^. 

Figure El presents the projection on Mx and x\ of the fitting result for the data in the 
whole region. The extracted yields for the signal components are N{^T~^^v) = 72 ± 11 
and N{p~l^v) = 59 ± 15, with the LCSR model used for the two signal PDF. 




Mx(GeV/c') Xb2(MJ XB2(Mp) 



FIG. 3: Projected Mx distribution(left) and distributions for the mass region of vr" {Mx < 
0.18GeV/c^, middle) and p~ (0.53 < Mx < l.lGeV/c^, right) in all region; open circle is data. 
Histogram components are 7r~£"'"z^ (red 135° hatch), p'^l'^u (blue 45° hatch), other Xu^-^v (green 
cross-hatch) and BB background (bluish green shaded). 

Figure m shows projections of the data, separated into three bins, 8 < < 16 GeV^/c^ 
and q^ > 16 GeV^/c^. Here the normalization of the other B^ X^i^u and the BB 
background components are fixed to those obtained in the above fitting for the whole q"^ 
region. Tableland Table ITTl summarize the extracted branching fractions with different FF- 
models for the 7r~i~^i' and p'i^u, respectively. The results are unfolded using an efficiency 
matrix that relates the true and reconstructed for the three intervals. We calculate the 
total branching fraction by taking sum of the branching fractions in the three q^ intervals. 

Table lllll summarizes the considered experimental systematic errors for the branching 
fractions. The total systematic error is the quadratic sum of all individual ones. A major 
contribution comes from the efficiency calibration with the B^ig D*^l^v sample, where 
8.3% originates from the statistics of the detected D*~C.^v decays and 4.9% from the error 
on B{B^ D*~i~^h') quoted in [i^. We consider the uncertainty in the number of B^B^ 
pairs; the ratio of B^B~ to B^B^ pairs (/+//o, 2.4%), the mixing parameter (x^, 1.0%) and 
the measured number of BB pairs {N^g, 0.5%). 
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FIG. 4: Projected Mx distribution in each q^ region; Open circle is data. Histogram components 
are tt~£~^u (red 135° hatch), p~£~^v (blue 45° hatch), other X^^.^i' (green cross-hatch) and BB 
background (bluish green shaded). 



TABLE I: Extracted branching fractions for B^ tt^I^u with different FF models in unit of 10~^; 
the total branching fraction and the partial branching fractions in three q'^ intervals, x^/''^ s-^id the 
probability of shows the quality of the fit of the FF shape to the extracted distribution. 



FF model 


Btotai 'B<8 Bs-16 13>i6 x^/n Prob. 


BaU'Ol 
ISGW II 
UKQCD 

Average 


1.77 ± 0.28 0.72 ± 0.15 0.57 ± 0.16 0.48 ± 0.17 2.4/2 30.3% 
1.75 ± 0.28 0.71 ± 0.15 0.60 ± 0.17 0.44 ± 0.16 1.9/2 37.8% 

1.75 ± 0.27 0.70 ± 0.15 0.60 ± 0.16 0.45 ± 0.16 0.4/2 81.0% 

1.76 ±0.28 0.71 ±0.15 0.59 ±0.16 0.46 ±0.17 - 



The systematic error due to the uncertainty on the inclusive branching fraction B{B — >■ 
Xuiu), which is used to constrain X~i~^i> background, is estimated by varying this parameter 
within ±1(7 of the error. The uncertainty on the BB background shape in our selection cut 
(iV^+ = 1, A^^o < 1) is studied in the simulation by randomly removing charged tracks and 
tt" according to the error in detection efficiency (1% for a charged track, 3% for vr''), and also 
by reassigning identified charged kaons as pions according to the uncertainty in the kaon 
identification efficiency (2%). The resultant change in the extracted branching fractions are 
assigned as a systematic error. We have seen significant change for p~C.^v. 

The FF model dependence of the extracted branching fractions has been studied by 
repeating the above fitting procedure by varying the FF model for the signal mode and also 
for the cross-feed mode {ttIv ^ piu). For the extracted B{B^ —>■ n'i'^i'), the FF model 
dependence is found to be small; the variation from the average is < 0.6% for the change of 
the iriu FF model and < 1.0% for the change of the piu FF mode. As for B{B^ —>■ p~t^v), 
significant dependence on the plv FF model is found; the variation is < 10.7% for the 
change of the plv FF model, while it is only < 1.3% for the change of the ttIv FF model. A 
similar tendency was observed in the CLEO analysis Q. At this stage, by adding linearly 
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TABLE II: Extracted branching fractions for B — > p^H.'^v with different FF models in unit of 10~ ; 
the total branching fraction and the partial branching fractions in three intervals, /n and the 
probability of shows the quality of the fit of the FF shape to the extracted distribution. 



FF model 


Btotai B<8 13>ie x^/n Prob. 


Ban'98 
ISGW II 
Melikhov 
UKQCD 

Average 


2.59 ± 0.73 0.24 ± 0.35 1.86 ± 0.51 0.50 ± 0.39 6.3/2 4.3% 
2.50 ± 0.86 0.18 ± 0.57 1.76 ± 0.48 0.56 ± 0.43 1.3/2 52.2% 
2.79 ± 0.83 0.20 ± 0.29 1.88 ± 0.51 0.72 ± 0.59 8.2/2 1.6% 
2.28 ± 0.65 0.22 ± 0.37 1.73 ± 0.46 0.33 ± 0.27 4.2/2 12.1% 
2.54 ±0.78 0.21 ±0.40 1.81 ±0.49 0.53 ±0.42 - 



TABLE III: Summary of systematic errors(%) of the branching fractions 



Source 






Tracking efficiency 


1 


1 


tt" reconstruction 




3 


Lepton identification 


2.1 


2.1 


Kaon identification 


2 


2 


D*ii' calibration 


9.8 


9.8 


Br{Xuiy) in the fitting 


0.2 


3.4 


BB background shape 


4.4 


31.5 


Nbb 


0.5 


0.5 


U/fo 


2.4 


2.4 


Xd 


1.0 


1.0 


total 


11.5 


33.5 



the maximum variations with the signal and cross-feed FF models, we conservatively assign 
1.6% (12%) as the systematic error due to the FF model dependence for the 7r~(p^)£"'"z/ 
total branching fraction. 

RESULTS 

The total branching fractions are obtained by taking the simple average of the values 
obtained from the FF models shown in Tables |l] and |n] 

B{B^ n-tu) = (1.76 ± 0.28 ± 0.20 ± 0.03) x 10-^ (2) 
B{B° p~tu) = (2.54 ± 0.78 ± 0.85 ± 0.30) x 10~^ (3) 

where the errors are statistical, experimental systematic, and systematic due to form-factor 
uncertainties. The obtained branching fractions are consistent with the existing measure- 
ments by CLEO for 5° ^ n'/p-i+u |6|] and BaBar for p~i+u 0, within the mea- 
surement uncertainties. 
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Figure IHl presents the obtained distributions for the two decay modes, overlaid with the 
best fits of FF shapes to the data. To be self-consistent, the shape of a particular FF model 
is fit to the distribution extracted with that FF model. The quality of the fit in terms of 

and the probability of , shown in Table H] and IHl may provide one way to discriminate 
among the models. At the present accuracy, we are unable to draw any conclusion on this 
point. 




FIG. 5: Extracted distrubution for the ^ 7r-£+zy(left) and 5° p-£+ i/(right) decays. Data 
points are shown for different FF models used to estimate the detection efficiency. Lines are for 
the best fit of the FF shapes to the obtained q^ distribution. 



We extract \Vub\ using the relation, 

\Vub\ 



where Tthy is the form-factor normalization, predicted from theories. In this paper, our 
major focus is on the \Vub\ determination based on the 7r~£'^z/ data and the form factor 
predicted by LQCD calculations. Since the current LQCD calculations are available only in 
the region > 16 GeV^/c^, we use the branching fraction in the high bin extracted with 
UKQCD; i3>i6 = (0.45 ± 0.16) x 10"! We use Tbo = 1.536 ± 0.014 ps for the 5° lifetime 

We apply tthy predicted by the FNAL [H, JLQCD (2J], APE as well as UKQCD 
calculations, as quoted by the CLEO analysis in 2003 0]. For the average of these results, 
the combined tthy = l-^^-oH 

± 0.47 ps ^ calculated by CLEO work is used. Here the 
errors are the statistical and the systematic in LQCD calculations, the latter including the 
quenching error of 15%. We obtain 

|Kfe|(;t>i6) = (3.90 ± 0.71 ± 0.23l°i2) ^ ^o-^ (5) 
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where the errors are statistical, experimental systematic, and theoretical. Figure IHl shows 
\Vub\ determined with these LQCD form factor predictions. Our result is compatible and 
consistent with that from CLEO in 2003. 

Recently, two preliminary results from the unquenched LQCD calculations were reported, 
FNAL'04 H and HPQCD [3]. Their reported fthy for the > 16 GeVVc^ region are 
Tthy = 1.96 ± 0.51 ± 0.39 ps^^ (FNAL'04) and 1.31 ± 0.13 ± 0.30 ps^^ (HPQCD). We also 
estimate \Vub\ using these calculations: 

|K6|Lt>i6) = (3.87 ±0.70 ±0.221°:^?) X 10-=^ (FNAL'04), (6) 



'«b|(g2>l6) 

|K6|(,t>i6) = (4.73 ± 0.85 ± 0.27l°:5^^) x 10"^ (HPQCD). (7) 



Finally, Figure [7| summarizes the Belle \Vub\ results from this work and from inclusive 
B — > Xuiiy measurements. The above reported results from 7T~i~^i' high data are 

in agreement with those from the inclusive B X„£z/ decays |^ 2^ 27 1. In the future, by 



accumulating more luminosity and with improvement in unquenched LQCD calculations, 
the present measurement will provide a more accurate determination of \Vub\ 



1 1 1 1 1 1 1 1 1 1 1 
Belle TTlv w/ LQCD (quenched) 


FNAL 


H • M 


JLQCD 


1— 1 • 1 1 


APE 


1 1 • 1 1 


UKQCD H 


• — ' 


Average 


Kl • 1-1 


(CLEO2003) 1 « 


Belle ttIv w/ LQCD (unquenched, preliminary) 


FNAL04 


HH # 1 1 


HPQCD 

1 1 1 1 1 


H • _ 

1 1 1 1 1 1 



0.002 0.004 0.006 



FIG. 6: \Vub\ determined using LQCD and the 7:~£~^i' branching fraction in the > 16GeV^/c^ 
region. For each data point, ztlo" range of the experimental errors (quadratic sum of statistical and 
systematic errors) and the range of the total error including the theoretical error are indicated. 
The CLEO result 6] with the averaged quenched LQCD formfactor is also shown for comparison. 
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\ 1 1 1 1 \ 1 ^ 1 1 r 

Belle Ttlv w/ LQCD (quenched) 



Quenched average 



Belle Tilv w/ LQCD (unquenched, preliminary) 
FNAL04 



HPQCD 



Belle X„lv 



(mj^,q ) w/ fuU-recon. 



(m^,q ) w/ v-recon. 



Pe end-point 



J_ 



J I L_ 



0.002 



0.004 



0.006 



FIG. 7: Summary of Belle \Vub\ determination. Results based on the exclusive — > n^i'^i' 
measurement (g^ > 16GeV^/c^) in this work with the averaged queched LQCD form factor (top) 
and with preliminary unqueched LQCD form factors (middle two) are shown, together with those 
based on our inclusive B — > Xuiu measurements (bottom three). The lower two inclusive results 
are rescaled using the recent determination of the fu factor j^l • 
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